We select cluster members using different membership cri- 
teria, to study the properties of pre-main sequence stars with 
or without circumstellar disks. We use intermediate resolution 
spectra including the Li 1 6707.8 A line to derive radial and rota- 
tional velocities, binarity and to measure the Equivalent Width 
of the lithium line; these results are combined with X-ray data 
to study the cluster membership. Optical-IR data and Ha spec- 
tra, these latter available for a subsample of our targets, are 
used to classify CTTS and WTTS and to compare the prop- 
erties of stars with and without disks. We find a total of 237 
certain members including 53 binaries. The rotational velocity 
distributions of stars with IR excesses are statistically differ- 
ent from that of stars without IR excesses, while the fraction 

■ of binaries with disks is significantly smaller than that of sin- 
(^) gle stars. Stars with evidence of accretion show circumstellar 
< ^ disks; youth of cluster members is confirmed by the lithium 

\ abundance consistent with the initial content. As indicated by 

■ the disk-locking picture, stars with disks in general have ro- 
/"■n tational velocities lower than stars without disks. Binaries in 

■ NGC 6530 seem have undergone significant disk evolution. 
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Abstract. Mechanisms regulating the evolution of pre-main sequence stars can be understood by studying stellar proper- 
ties such as rotation, disk accretion, internal mixing and binarity. To investigate such properties, we studied a sample of 332 
candidate members of the massive and populous star forming region NGC 6530. 



1. Introduction 

Very young open clusters are crucial systems for studying star 
formation mechanisms in different conditions (density, metal- 
licity, age). The present study is part of an ongoing project to 
investigate the properties of NGC 6530, a star formation re- 
gion located in front of the giant molecular cloud M8 (Lagoon 
Nebula) and about 125 pc from the Sun . With stars having a 
median age of 2.3 Myr JPrisinzano et all 120051) . NGC 6530 is 
very suitable for investigating the stellar and disk evolution of 
pre-main sequence (PMS) stars. 

Several papers have bee n devote d to this clu ster (Walker, 
1957| iThel Il96f)t Ivan den Ancker eTail Il997t ISun g et all 
2000: lRauw et all l20021: ts low-mass population has b een re - 
cently investigated by using X-ray data (Damia ni et all 12 004). 
deep BVI photometry (d own to V ~ 23 , Prisinza no et all 2005) 
and 2MASS IR data JDamiani etail 120061) . In these recent 
studies, more than 1100 candidate cluster members have been 
identified using X-ray data and/or optical-IR photometry. 

One of the most challenging tasks in studying open clus- 
ters is to identify certain members, especially in case of PMS 
stars which are located in a broad region of the optical color- 
magnitude diagram. Spectra around the lithium and Ha lines, 
such as those analyzed in this work, are very useful to derive 
different and independent membership criteria, based on the ra- 
dial velocities and the lithium line. 

In addition, spectroscopic observations are useful to iden- 
tify binaries, to measure projected rotational velocities, lithium 
equivalent widths and to distinguish stars with accretion disk 
phenomena by using the Ha line. These are crucial pieces of 
information to study the properties of PMS stars and the evolu- 
tion of their disks. 

PMS stars with accretion disks in which the material falls 
onto the stellar surface are usually classified as CTTS. These 

Send offprint requests to: e-mail: loredanaSastropa . inaf . it 



objects can be spectroscopically identified since the heating of 
the in-falling gas produces broad emission lines, such as the 
Ha line. On the contrary, stars without accretion phenomena, 
which show narrow chromospheric emission lines, are classi- 
fied as WTTS. 

In the case of a lack of spectroscopic observations indicat- 
ing accretion phenomena, stars with a disk, accreting or not, 
can be identified using their optical-IR photometric proper- 
ties, as already discussed in several papers (see for example 
ISicilia-Aguilar et all l2005t [Damiani et all I2006T) . circumstel- 
lar disks produce IR excesses in the spectral energy distribution 
(SED) with respect to ordinary stars. Stars with IR excesses, 
in a small region around the Hourglass ne bula in the M8 star - 
forming region, have also been identified bv l Arias et all ll2006); 
these stars are classified as low- and intermediate-mass pre- 
main-sequence star candidates. 

The excesses can be quantified using a combination of 
optical-IR colors, which become useful means to distinguish 
stars with or without a disk. However, while the selection of 
CTTS by means of the Ha line is in general agreement with 
the s election of stars with disk de tected by optical-IR photom- 
etry JSicilia-Aguilar etall Eo05). a clear distinction between 
CTTS and WTTS cannot be obtained using only photometric 
data, since not all stars with a disk show accretion phenomena. 

Circumstellar disks in PMS stars seem to be strictly 
related to the angular momentum transf er mechan is ms, a s 
recently found, for ex a mple, bv j He rbst & Mundtl §005); 
iDullemond et all J2006I) : iRebull et all J2006I) . Therefore, the 
comparison between the rotational velocities of single PMS 
stars surrounded by a disk and those of PMS stars which lack 
a disk, allows us understanding of how the angular momentum 
can be regulated by circumstellar disks. Using a statistically 
large sample of certain single cluster members and robust IR 
excess indicators is crucial to test if stars surrounded by disk 
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are slower than stars wit hout disk, as asserted in the disk lock- 
ing model jKoenig]Lll99ll) . 

It is widely believed that the majority of stars form in bi- 
nary systems, but the mechanisms of the evolution of their cir- 
cumstellar disks and their contributio n to the IR fluxe s from 
the companion are not clear JGras-Velazauez & Ravi l2005t 
iMonin et all I2006I) . Evidence of how disk properties vary in 
the binary environment and with binary parameters, such as 
separation and mass, has not been found yet, since high resolu- 
tion photometry and spectroscopy have only become available 
in the last decade. 

Another controversial question is whether the star forma- 
tion process is instantaneous or not. Since the stars deplete their 
initial lithium content in the first few million years, lithium 
abundances in young stars can be used as age indicators to 
study the ag e spread and hence the duration of the star forma- 
tion process JPalla et all 120051) . However, accurate derivations 
of lithium abundances are hampered by several difficulties re- 
lated to the effective temperature measures and to the theoret- 
ical models, especially for stars cooler than 4000 K, where the 
physics of the stellar atmosphere is not well known. 

A combined study of X-ray data, optical-IR photometry 
and spectroscopic data is presented in this paper, with the aim 
to investigate the main properties of a sample of young stars 
in NGC 6530, such as relations between disk presence and an- 
gular momentum, accretion phenomena, binarity and lithium 
abundance. In Section 2, we present the criteria adopted to se- 
lect our targets and the spectroscopic data; in Section 3 we de- 
scribe the data analysis and in particular the technique used to 
derive radial and rotational velocities and to distinguish single 
stars from binaries; the methods adopted to measure the lithium 
Equivalent Width and to classify CTTS and WTTS stars based 
on the He line are also described. In Section 4, we present the 
analysis of stars with or without optical-IR excesses, the veil- 
ing correction and the adopted membership criteria. Properties 
related to the rotational velocities, accretion phenomena and 
binarity of IR excess stars are presented in Section 5; the dis- 
cussion on the lithium abundance is presented in Section 6 and 
final results are summarized in Section 7. 

2. Target selection and observations 

A list of candidate members of the star forming region 
NGC 6530 suitable for spectroscopic investigations was ex- 
tracted from the optical catalog obtained from images taken 
with the WFI ca mera of the ESQ 2.2m Telescope and has 
been presented by Prisi nzano et alJ d2005h . As already in that 
latte r paper, most o f the stars with an X-ray counterpart 
( Da miani et alll2004l) are located in a quite well defined region 
of the V vs. V-I color-magnitude diagram (CMD), since most 
of them are cluster members. This region of the CMD identi- 
fies the PMS cluster locus. We selected our targets taking stars 
falling in this region of the CMD, with V magnitude between 
14.0 and 18.2. Most of the selected targets (about 72%) are X- 
ray active stars, while the remaining ones are objects with sim- 
ilar photometric properties. Fig.^ shows the CMD of our tar- 
gets where X-ray detected objects are indicated with large bul- 
lets. The solar metallicity ZAMS (dashed line) and isochrones 
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Fig. 1 V vs. V-I color-magnitude diagram of the stars in the 
NGC 6530 field observed with Flames. X-ray detected targets 
are drawn as large bullets. The dashed and solid lines are, re- 
spectively, the solar metallicity ZAM S and isoc hrones of 1 .0, 
3.0 and 10.0 Myr, computed bv lSiess et all feOOOh . 

(solid lines) of 1.0, 3.0 and 10.0 Myr, computed bvlSiess et all 
( 2000]) are also drawn, using the iKenvon & Hartmanr]7l995F 
transformation s, the r eddening value of E(B - V) = 0.35 given 
bv ISung et al] J2 000l) and the cluster distance d=1250pc de- 
rived in lPrisinzano et alJ d2005l) . 

A list of our selected targets is presented in Tabled where 
columns 1 to 6 give the celestial coordinates, col. 7 g ives 
the spectrum number, col. 8 is the Damiani et al. (200 41 X- 
ray identification number, col. 9 is the Prisinzano - e^aljj|2005j) 
WFI id entification number, col.s 10-15 are the Prisin zano et alJ 
(2005) BVI magnitudes and their photometric uncertainties, 
col.s 16-21 are the IR magnit udes and photome tric uncertain- 
ties 1 from the 2MASS catalog dCutri et alll2003l) while col.s 22 
and 2 3 give ages and masses as estimated in Prisinzano et al. 
(2005$, from theoretical isochrones and tracks dSiess et all 
120001) 

The targets were observed with the GIRAFFE-FLAMES 
multi-object fiber spectrograph of the ESO-VLT Kueyen 
Telescope (UT2). We observed 332 objects with the Medusa 
fibre system in combination with the 316 lines/mm grating 
and the setup HR15 which gives a resolution R=19300 in the 
wavelength range [6606, 6905] A, comprising the Li 1 6707.8 A 
line. Furthermore, 115 of our targets were also observed us- 
ing the setup HR14, which gives a resolution R=28800 in the 
wavelength range [6383, 6626] A, comprising the Ha line. 

1 'NaN' indicates magnitudes rejected since labeled as 95% con- 
fidence upper limits i.e. sources not detected or inconsistently de- 
blended in that band. 



Table 1 List of stars observed with FLAMES. Sp is the spectrum number, IDx is the source identification number of Da miani et alJ (2006): IDwfi is the Prisinza no et al] (12005 > 
identification number; B, V, I and eB, eV and el are the magnitudes and their photometric uncertainties presented in IPrisinzanoetal.1 J2005h . while JHK and eJ, eH and eK 
are the 2MASS magnitudes and their photometric uncertainties; "NaN" indicate magnitude s rejected since labeled as 95% confidence upper limits i.e. sources not detected or 
inconsistently deblended in that band. Ages and masses are those estimated in Prisinza no et al.l (12005). The complete table is available in electronic format at CDS. 



RA(2000) 


Dec(2000) 


Sp ID X IDwfi 


B 


eB 


V 


eV 


J 


el 


J 


eJ 


H 


eH 


K 


eK Age 


Mass 


h m s 


dm s 


























Myr 


M e 



18 4 


24.59 


-24 


19 


24.69 


2002 


493 


25079 


17.280 


0.002 


15.836 


0.002 


13.985 


0.002 


12.414 


0.028 


11.447 


0.028 


10.888 


0.028 


0.77 


1.05 


18 4 


24.62 


-24 


20 


4.93 


2003 


495 


23868 


17.118 


0.004 


15.767 


0.001 


14.082 


0.002 


12.832 


0.032 


12.067 


0.032 


11.833 


0.032 


1.35 


1.38 


18 4 


7.60 


-24 


12 


21.43 


2004 




39006 


19.048 


0.010 


17.455 


0.005 


15.317 


0.005 


NaN 


NaN 


NaN 


NaN 


12.562 


NaN 


1.62 


0.67 


18 4 


23.67 


-24 


17 


29.10 


2005 


470 


27805 


18.486 


0.003 


16.822 


0.003 


14.573 


0.003 


12.948 


0.050 


12.090 


0.050 


11.789 


0.050 


0.62 


0.56 


18 4 


16.11 


-24 


19 


52.26 


2006 


301 


24253 


17.227 


0.002 


15.790 


0.002 


13.893 


0.005 


12.373 


0.036 


11.449 


0.036 


10.829 


0.036 


0.64 


0.98 


18 3 


59.44 


-24 


12 


23.42 


2008 


70 


38914 


16.939 


0.004 


15.708 


0.003 


14.146 


0.004 


13.083 


0.047 


NaN 


NaN 


NaN 


NaN 


2.31 


1.66 


18 4 


13.38 


-24 


15 


13.49 


2010 


243 


32646 


18.545 


0.004 


16.995 


0.005 


14.905 


0.002 


13.366 


0.058 


NaN 


NaN 


12.246 


0.058 


1.27 


0.70 


18 4 


15.79 


-24 


21 


27.34 


2012 


295 


21339 


17.666 


0.002 


16.335 


0.002 


14.610 


0.002 


13.314 


0.026 


12.608 


0.026 


12.350 


0.026 


2.35 


1.29 


18 3 


57.28 


-24 


16 


9.83 


2013 


49 


30517 


17.468 


0.003 


16.317 


0.002 


14.871 


0.003 


13.507 


0.072 


12.676 


0.072 


12.136 


0.072 


9.27 


1.38 


18 4 


21.79 


-24 


14 


3.91 


2014 


426 


35254 


15.237 


0.004 


14.229 


0.002 


13.033 


0.007 


12.117 


0.022 


11.644 


0.022 


11.486 


0.022 


6.97 


1.83 
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These observations were made on May 27, 2003 and are part 
of the Guaranteed Time Observations of the Ital-FLAMES 
Consortium. All targets were observed using three different 
FLAMES fiber configurations, as described in the Log-book 
in Table|2] Spectra taken by using the configuration ID 090229 
were summed. 

The data were reduced using the GIRAFFE Base-Line 
Data Reduction Software (Version 1 .08) developed at Geneva 
Observatory and relative calibration files (Version 2.0). 

Fig-El snows two examples of reduced spectra around the 
Li I 6707.8 A line for an X-ray detected star {solid line) and 
for an X-ray undetected star {dashed line). Both spectra show 
several Fe I an d Al I lines , the n ebular sulphur emission line 
at X = 6717 A jLada et allll976l) and the Ca I 6717.69 A ab- 
sorption line; only the upper spectrum shows a very strong Li 
feature. 

The sky contribution from the spectra obtained around the 
Li 1 6 707.8 A line was subtracted by using a median, computed 
as in I.Teffries & Oliveiral ll2005l) . of the 15 sky fiber spectra 
available for each configuration and obtained during the same 
night. The sky contribution is mainly dominated by emission 
of the H II region, which is spatially varying as suggested from 
the strong nebular emission lines which can have intensity vari- 
ations up to 90%. For this reason, the strong nebular emission 
lines have not been successfully subtracted, while a good sub- 
traction of the sky continuum was achieved in the vicinity of 
lines that will be analyzed in this work. 

Due to the strong and variable contribution of the nebular 
emission Ha line, we cannot subtract the sky level from our 
spectra in the Ha region. The analysis of these spectra will be 
described in Sect. 13.31 



Table 2 Log of the observations. 



Config. 


Exp. Time 


Grating 


# stars 


spectrum number 


V 


ID 


[sec] 






range 


range 


214312 


5400.0 


HR15 


107 


[1002-1135] 


[14.0-18.2] 


174741 


2808.9 


HR15 


114 


[2002-2135] 


[14.0-17.7] 


174741 


1211.0 


HR14 


114 


[2002-2135] 


[14.0-17.7] 


090229 


2697.1 


HR15 


111 


[3002-3135] 


[14.0-18.2] 


090229 


3000.0 


HR15 


111 


[3002-3135] 


[14.0-18.2] 



3. Data analysis 

3. 1. Radial and rotational velocities and binarity 

We computed radial velocities (RV) and rotational velocities 
of our targets using the IRAF task FXCOR, which performs 
a Fourier Cross-Correlati on of each spect r um w ith a template 
spectrum, as described in iTonrv & Davisl Jl979l) . We used as 
template the spectrum of a relatively bright star {V - 14.5) of 
our sample (Sp 1083), which is regarded as a reliable cluster 
member since it is a photometric candidate PMS star, with X- 
ray emission and a very strong lithium line {Fig$2]solid line), 
which is a feature of very young stars. Among the stars with 
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6680 6690 6700 6710 6720 6730 

X [A] 

Fig. 2 Examples of spectra around the Li I 6707.8 A line for an 
X-ray detected star, Sp 1083 {solid line) and for an X-ray un- 
detected star, Sp 1052 {dashed line). Both spectra show several 
Fe I and Al I lines, the nebular sulphur emission feature and the 
Ca I 6717 A absorption line; only the upper spectrum shows a 
very strong Li feature. 

these characteristics, this star also shows a spectrum with the 
narrowest lines, having a FWHM of about 20km/s, which is 
slightly larger than the instrumental spectral resolution of our 
spectra of about 17km/s, as estimated by our lamp spectra; this 
latter value sets a lower limit to the rotational velocities we can 
determine. 

For each star of our sample, we cross-correlated the spec- 
trum with the template over the wavelength range between 
about 6610 and 6820 A, which contains no significant telluric 
absorption lines. Within this wavelength range we selected only 
regions without sky or nebula emission lines. 

The position of the cross-correlation peak is a measure of 
the relative RV of each star with respect to the RV of the tem- 
plate star, since it is a measure of the Doppler shift between the 
two cross-correlated spectra. 

For 53 stars of our sample, the peak of the cross-correlation 
function is not symmetric and can be fitted with two over- 
lapping Gaussians rather than with a single one, as shown in 
the two examples of Fig. [3] where the relative radial velocities 
of the two components is also indicated. In these cases, the 
cross-correlation function is indicative of a double-lined spec- 
troscopic binary (SB2). For these objects, however, we cannot 
measure the RV of the binary systems, since our spectra were 
obtained during a single night. 

To obtain the RV in the heliocentric system we proceeded 
as follows: a) we used the IRAF task RVIDLINES to derive 
the observed RV of the star used as template, i.e. the wave- 
length shift in spectral lines relative to rest wavelengths; we 
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Fig. 3 Example of peaks of the cross-correlation functions pro- 
duced by cross-correlating two binaries against our template 
star. Relative radial velocities of the two components are indi- 
cated in each panel. 



find that the observed RV of this star is -11.85 + 0.56 km/s, 
based on 48 lines; b) using the log-book of the observations 
and the IRAF task RVCORRECT, we find that the heliocentric 
correction to the observed RV is equal to 12.310 + 0.005 km/s; 
c) we computed the heliocentric RV of the template star, which 
is equal to 0.47 + 0.56 km/s, and from this, the heliocentric RV 
of all stars (see Tablel^- 
The histogram in Fig. is the RV distribution of the whole 
sample of 279 spectroscopically observed single stars; the his- 
togram shows a significant peak indicating the presence of 
the cluster. In order to distinguish cluster members from field 
stars, we fitted this distribution with a double Gaussian us- 
ing the "maximum likelihood fitting". We find that the cluster 
Gaussian is centered on 0.5+0.3 km/s with a standard deviation 
cr =3.2+0.2 km/s {solid line in Fig.0}, while the much broader 
field star RV distribution shows apeak at -45.2+ 10.5 km/s, with 
a standard deviation of cr =101.0±7.5km/s (dashed line in 
Fig-Ell- These values are affected by large uncertainties due to 
the poor statistics of the field star RV distribution, reflecting 



our bias towards likely cluster members, selected from their 
position in the CMD and the X-ray flux. 

The total number of possible cluster members within +3cr 
of the cluster RV distribution is 190 including ~6 contaminat- 
ing field stars, as computed from the Gaussian distribution of 
the field stars. 

Fig-0shows the comparison of the RV of the whole sample 
of single observed stars (279), with the RV distribution of the 
187 X-ray detected single stars and with the RV distribution of 
the 92 X-ray undetected stars. The latter shows a weak peak at 
the cluster RV which is due to the X-ray selection incomplete- 
ness, as will be discussed in Section l4~3l 

L 1 1 1 1 I I I _ 



80 - 



40 - 



20 - 




-200 -100 100 

Hel. Rad. Vel. [km/s] 

Fig. 4 RV distribution (solid histogram) of the whole sample 
of single stars which includes 279 objects. The Gaussian RV 
distributions for the cluster (solid line) and field (dashed line), 
obtained from the double Gaussian fitting, are also shown. 

The cross-correlation function allowed us to measure also 
the projected rotational velocities v sin(i) of our targets, since 
the width of the cross-correlation peak is a function of the ro- 
tational broadening of the spectrum lines. As discussed before, 
our template is the spectrum of a slow rotation star and then the 
width of the cross-correlation peak with each target is a mea- 
sure of the rotational velocity of each target with respect to the 
rotational velocity of the template. 

Fig-ED shows 4 examples of spectra with different broad- 
ening of the absorption lines and therefore different rotational 
velocities. The upper spectrum is that of the slow rotation tem- 
plate star. 

To calibrate the relationship between the width of the cross- 
correlation peak and the v sin(i) values, we computed sy nthetic 
spectra at different v sin(i), using the ATLAS9 lKuruc3 Jl993l) 
model atmospheres, assuming the GIRAFFE instrumental res- 
olution R= 19000, temperature T e ft=4600 K (estimated from 
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Fig. 5 Zoom of RV distribution of the whole sample of the 279 
single stars, compared with the RV distribution of the 187 sin- 
gle X-ray detected stars and with the RV distribution of the 92 
single X-ray undetected stars. 

V - I, see Sect. l3.2t . log g = 4.0, microturbulence £ = 1 km/s, 
macroturbulence 2 km/s and solar metallicity, which are the ex- 
pected parameters of the template star. 

The obtained synthetic spectra were cross-correlated with 
our template and the FWHM of the cross-correlation peak 
was measured. To check the sensitivity of this relationship on 
the temperature, two additional sets of synthetic spectra were 
computed at T e ff=4000K and T e ff=5000K, which are typical 
temperature values for our targets, estimated as described in 
Sectionl3~2l 

The resulting relations between the FWHM of the cross- 
correlation peak and the corresponding vsin(i) for the three 
considered temperatures are shown in Fig. [7] We note that the 
three relations are not temperature dependent up to v sin(i)^ 
60 km/s; for vsin(i) > 60 km/s, the discrepancies are, however, 
consistent within the typical errors inferred for the rotational 
velocities measured for our targets (see below in this section). 
The lower limit of the v sin(i) values we can measure is ap- 
proximately set by the instrumental resolution, equal to about 
17 km/s, which corresponds to the flattening of the above rela- 
tions of calibrations. 

Using the relationship corresponding to T e ff=4600K, 
which is the expected temperature of our template, we interpo- 
lated the FWHM of the cross-correlation peak with the derived 
calibration curve to estimate the v sin(i) values. Rotational ve- 
locities were derived only for single stars, since the double 
peak of the cross-correlation function in SB2 stars cannot be 
accurately separated into the two componen ts to obtain reli- 
able vsin(i) values. As in lRhode et alJd200lh . the errors of the 
rotational velocities were taken equal to + v ""^' ) , where r is the 






6680 



Sp 1079 v sin(i)= 24.22± 1.32 



6700 6720 
A [A] 

Sp 1024 v sin(i)= 62.84± 6.46 



6700 6720 
A [A] 

Sp 2055 v sin(i) = 145.11 ±13.73 



6700 6720 
A [A] 



Fig. 6 Examples of spectra with different broadening of the 
spectrum lines and therefore with different rotational veloci- 
ties. The upper spectrum is that of the slow rotation template 
star. 



iTonrv & Davisl ( ll97 9) parameter, which is a measure of the sig- 
nal to noise ratio of the cross-correlation peak. 

Of the 279 observed single stars, 140 have v sin(i) less than 
or equal to our estimated actual limit (~ 17 km/s) for the rota- 
tional velocities. For 7 stars we estimated vsin(i) > 140 km/s, 
but the signal to noise ratio of the cross-correlation peak for 
these stars is very low (r < 10) and therefore it is not possi- 
ble to distinguish whether they are binaries and/or fast rotators. 
The large v sin(i) values obtained for these stars are therefore to 
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Table 3 Results of the spectral analysis obtained in this paper for the star listed in Tabled See Sections l3~2l 13 . 31 |4~T1 and I4.3l for 
explanations. The complete table is available in electronic format at CDS. 
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fc C=CTTS, C?= possible CTTS and W=WTTS 

c PT=Hci' profile type {see text) 




Fig. 7 Calibration curve used for the vsin(i) estimation from 
the width of the cross-correlation peak. 



Fig. 8 vsin(i) values vs. stellar masses for single stars. Cluster 
candidate members based on their X-ray detection and/or RV 
values are indicated. 



be taken with caution. Fig.|S]shows the v sin(i) values as a func- 
tion of the stellar mass, for the X-ray detected stars {bullets) and 
for the RV candidate members {squares). The results obtained 
as described in this section are given in Table [3] where col. 1 
gives the spectrum number, col. 2 indicates whether the star is 
SB2 (label "Y") or possible SB2 (label "?") while col.s 3-6 
give the heliocentric RV values and the v sin(i) values, with the 
estimated errors. 



3.2. Lithium equivalent widths 

To measure the equivalent width (EW) of the lithium line of 
our targets, we first performed the continuum normalization in 
the region around the Li I 6707.8 A line. We used the IRAF 
task CONTINUUM, considering a small line-free spectral range 
(< 10 A) around the Li 1 6707.8 A line, a second order Legendre 
function and a variable residual rejection limit. Due to the stel- 
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lar rotation, the FWHM of the line can vary from about 0.3 A 
up to 2.8 A. For this reason, both the line-free spectral range 
and the rejection limit were interactively chosen for each case, 
based on the visual inspection of the fitting results. 

Equivalent widths of the Li I 6707.8 A line were measured 
with the IRAF task SPLOT; since the EW of the lithium line 
is in most cases greater than 100 mA, we assumed a Voigt pro- 
file for the majority of our targets which show a very strong 
lithium line; a Gaussian profile was considered only in the few 
cases where the EW of the line was smaller than about 100 A. 
The EWs were not corrected for the contribution due to the 
nearby Fe I line at 6707.44 A blended with the Li line, since 
for our targets this contribution is negligibl e with respect to 
the lit hium. Following the method given in Soder blom et alJ 
l l 1993b . for the stars of our sample the EW of the Fe can range 
from 10 to 20 mA, which implies a systematic contribution of a 
few percent for our EW measures, that are typically larger than 
200 mA. 

In order to estimate uncertainties in the EW measures, we 
repeated four times both the continuum normalization and the 
EW measures. The maximum error distribution of such mea- 
sures does not depend on the S/N nor on the FWHM, but 
slightly increases with the mean value, up to a value of about 
0.15 A for EW~ 0.6 A. We then considered as final EW values 
the mean of the four measures and as typical error the stan- 
dard deviation of the maximum errors of all measures, which 
is 0.035 A for all EW. 

Values of the lithium EW for SB2 stars were computed with 
respect to the combined continuum which is the sum of the 
continuum contributions of both components. Lacking infor- 
mation on the individual spectral types of the two components, 
we cannot separate their continuum contributions and there- 
fore we cannot obtain individual estimates of the lithium EW 
for these stars. 

The EW of the lithium line in young stars depends strongly 
on the star temperature. Lacking information of spectral types, 
accurate temperatures of our targets are not available. For this 
reason, temperature estimates have been derived from photom- 
etry as follows: we dereddened observed (V - 1) colors by sub- 
tracting the mean re ddening E(V - I) = 0.46, obtained from 
the E(B - V) = 0.35 JSung et allkOOOl) and the reddening law 
E(V - 7) = l.3 E(B -V) = 0.35 dMunari & Carrarol fl996). As 
already done in Prisinza no et al.l (12005). the dereddened (V—I)o 
colors w ere converted into effective t emperatures using the re- 
lation of Kenvon & Hartmann lll995l) . 

The errors on effective temperatures strongly depend on the 
unknown individual reddening, which is the largest source of 
error. To estimate temperature errors, we considered four typ- 
ical observed V - I values (1.25, 1.70, 2.15, 2.80) which we 
dereddened using the redde ning range E(B - V) = [0.25, 0.50], 
found bv lSungetalJ lfeoOO). Typical temperature errors are in- 
dicated as horizontal bars in Fig.|9j where the lithium EWs of 
all single members are shown as a function of temperature. 

In order to verify the young age of our targets, the lithium 
EWs of older stars in the Pleiade s cluster (100 Myr) were also 
plotted (Soderblo m et all Il993l) . The comparison of the two 
samples of stars clearly shows that the EW of most of our tar- 
gets is significantly higher than those of stars of similar temper- 



ature in the Pleiades. This is consistent with theoretical mod- 
els that predict a co smic abundance of lithiu m in s t ars younger 
than a bout 10 Myr fo'Antona & M azzitelhl ll997tlsiess et all 
2000). 
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Fig. 9 Lithium EW for candidate single members of NGC 6530. 
Lithium EW of the older Pleiades stars are also shown for com- 
parison. Horizontal bars indicate expected temperature errors, 
while the dashed line indicates the limit adopted to define can- 
didate Li depleted stars. The 8 stars with evidence of Li deple- 
tion and/or veiling effects are also indicated. 

EW estimates of the lithium line for our targets are reported 
in Table|3] where col. 7 gives the temperature of the stars, col. 8 
gives the EW of the lithium for single stars or for the first com- 
ponent, if the object is identified as SB2; in these latter cases, 
the lithium EW of the second component is given in col. 9. The 
label "NaN" indicates that either the line is not present or it 
cannot be measured. 

3.3. Ha observations 

All our targets are projected onto the Lagoon Nebula, a very 
bright HII region, which produces a strong and spatially non 
uniform emission line background. This means that all spec- 
tra of our targets include a contribution of a spatially variable 
Ha emission, due to the Nebula itself. Since this contribution 
cannot be individually measured and subtracted from the stel- 
lar spectra, we are not able to measure the purely stellar EW of 
the Ha line. 

By inspection of the Ha spectral region, however, it is obvi- 
ous that a number of stars show broad lines, while the Ha from 
the Nebula is very narrow. This allows us to identify these lat- 
ter stars as CTTS, since their Ha line indicates gas motion in 
extended envelopes. 
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Fig. 10 Examples of Ha line profile as a function of the 
Doppler velocity. Note that fluxes are shown on a log scale. 



To select stars with signatures of accretion we measured 
the Full Width at Zero Intensity (FWZI) of the Ha line for 
the subsample of 115 stars for which the spectrum around the 



Ha region is available. The measures are reported in Col. 10 of 
Tabled the errors depend on the signal to noise ratio and are 
of the order of a few A. 

Different kinds of Ha line profi les are found a nd classified 
following the notation described in Bertout ( 1984) (see col. 19 
of Table[3j; some examples are reported in Fig.^H As summa- 
rized in the tree diagram shown in Fig.[^ of all the 95 sin- 
gle stars for which the Ha spectrum is available, 71 are cer- 
tain members while 2 are possible members. Among the cer- 
tain members, 3 1 show signatures of accretion and we classify 
them as CTTS; in particular we have: 1 star (Sp 2106) show- 
ing a type I P Cygni profile (label 'PCI'); 28 stars showing 
from one to three broadened emission components, as for ex- 
ample Sp 2002 and 2060 of these objects are classi- 
fied as type III P Cygni profile (label 'PCIII'); finally, 2 stars 
have inverse P Cygni profiles, as Sp 2042 (label 'YY Ori'). In 
addition, we have 9 stars with slight broad wings which can 
be due to "normal" chromospheric emission or to low accre- 
tion; we label these objects as "CE?" and we classify them as 
possible CTTS or "CTTS?". On the other hand, we have a to- 
tal of 3 1 member stars which do not show evidence of accre- 
tion: 6 of them show a broad absorption Ha line, with the core 
filled of the narrow nebular emission Ha, that we label as ' AE' 
(Sp 2076 is an example), while the remaining 25 stars show 
a quite symmetric and "narrow" line profile, likely originating 
only from the nebular emission; we label these objects as 'EG' 
and an example is the star Sp 2135. We classify these 31 ob- 
jects as WTTS. The 2 possible members show a Ha line profile 
labeled as 'AE'. 

These results are also reported in Table[5] where col. 18 in- 
dicates if they are classified as CTTS (label 'C') or possible 
CTTS (label 'C?') or WTTS (label 'W'); col. 19 indicates the 
profile type of the Ha line. 

4. Analysis 

4. 1. Stars with infrared excesses 

IR excesses are indicators of presence of circumstellar disks, 
typical of T Tauri stars . A study of IR-excess stars in NGC 6530 
dDamiani et all 2006) has shown that different optic al-IR ex- 
cess indices are related to different kinds of spectral energy 
distributions. 

We use, here, the extinction-free indices 

gjHHK = (J - H) - E(J - H)/E(H - K) x (H - K), 

gviJK = (J - K) - E(J - K)/E(V - I) x (V - I), (1) 

Qvnj = (I - J) - E(I - J)/E(V - 1) x (V - 1), 

Gbvjk = (J - K) - E(J - K)/E(B - V) x (B - V), 

already defined in Da miani et alJ d2006l) . to select, in our sam- 
ple, stars with IR excesses in at least one of the four indices. 
The iRieke & Lebofskvl dl985l) interstellar extinction law has 
been used for the JHK magnitudes. We consider only objects 
with an IR magnitude error smaller than 0. 1 mag. To avoid non- 
member contaminating stars we consider only objects that are 
members according to at least two of our membership criteria, 
based on X-ray detection, RV and presence of the lithium line 
(see Section l4~3l . 
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Fig. 1 1 Tree diagram showing the classification of the 95 single stars with different kinds of Ha profiles. Labels are described in 
Sect-IP 

95 singles 
w/ Ho- 



1 PCI 




22 NM 



31 WTTS 
6 AE 25 EG 



28 PCIII 2 inv PC 



Fig.^]shows our computed optical-IR indices as a function 
of one optical or IR color, for our candidate single members; 
different symbols indicate the membership criterion. Typical 
propagated errors are plotted in e ach panel, w here the so lid lin e 
indicates normal star colors from Ken von & Hartmam] |fl995). 
for Eb-v - 0.5, which is the maximum reddening estimated for 
the cluster. 

We define "IR excess stars", the objects falling in the gray 
shade regions of Fig.^1 i.e . stars with indic es exceeding the 
fiducial limits defined in Da miani et al.l (120061) and indicated in 
each panel by the dashed line. The selection has been made 
regardless of photometric error. The stars that we have classi- 
fied as certain members (see Section 1431 having IR excesses 
with l-cr error of the index compatible with normal colors, are 
13 from a total of 81. Nevertheless, only a fraction of these 13 
stars could not be stars with IR excesses. 

The number of stars with each of our membership indica- 
tors that are also disk candidates based on each of the IR excess 
indices is given in Table0] For each membership indicator (X, 
RV, Li), the last two lines of Tab le|4] give also the number of 
stars showing IR excess based on at least one of the four in- 
dices (Q xxxx ) and the total number of stars with the given mem- 
bership indicator. 

Stars with and without IR excesses are reported in Table[3] 
where col.s 1 1-13 indicate whether the star has an X-ray coun- 
terpart and/or if it is member based on the RV and/or if the 
lithium line is present while col.s 14-17 indicate whether the 
star shows IR excess in the corresponding Q index. Properties 
of the stars with or without IR excesses are discussed in 
Section|5] 

4.2. Veiling 

As in |Pallaetal] lEo05). we estimated the amount of spectral 
veiling affecting our stars, using the three strong lines of V I 
at 6624.8 A, Nil at 6643.6 A and Fel at 6663.4 A. We cannot 
compare the EW of these lines in our targets with those mea- 
sured in stars of IC2391 and IC 2602 jRandichet all 1200 lh 
and used by Palla et al. (2005), since the spectral resolution 
is R ~ 43 800, significantly higher than that of our spectra 



Table 4 Number of stars with each of our member indicators 
(X-ray detection, RV measurements and EW of the Lithium 
line), which are also disk candidates based on the IR excess 
indices. For each membership indicator (X, RV, Li), the last two 
lines give also the number of stars showing IR excess based on 
at least one of the four indices (Q XX xx) and the total number of 
stars with the given membership indicator. 





X 


RV 


Li 


QjHHK 


26 


32 


29 


QviJK 


67 


76 


74 


QviIJ 


12 


14 


13 


Qbvjk 


68 


76 


74 


Qxxxx 


76 


84 


81 


Tot. 


187 


190 


195 



(R ~ 19 300). To give a consistent estimate of veiling in our 
targets, we compared the EW of the selected lines with those 
of WTTS candidates selected in our sample. WTTS candidates 
were chosen as those objects that are X-ray detected, members 
based on the RV, showing the lithium line, and without IR ex- 
cesses with respect to the Q reddening-free indices defined in 
Section ETTI This condition ensures that we select stars without 
strong disks and likely not affected by veiling. To avoid large 
errors due to rotational broadening in the EW measures of the 
three lines, we considered only slow rotation members (v sin(i) 
< 30km/s). 

Panels a, b and c of Fig. ^] show the EWs of the three lines 
measured for all our targets that are members according to at 
least one criterion. Our WTTS candidates are also indicated. 
We do not have a significant number of WTTS for temperatures 
higher than ~4700 K and therefore we can estimate the veiling 
contribution only for stars cooler than this temperature. 

For each of the three lines and each candidate, we com- 
puted the veiling correction as r = EWwTTs/EW st ars - 1 where 
EWwtts has been obtained from a linear fitting of the EW of 
the WTTS, for the V I line, while it has been taken as the me- 
dian value of the EW of the WTTS, for the Ni I line and the Fe I 
line, for which no trend of the EW as a function of temperature 
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Fig. 12 Extinction-free indices defined in Damia ni et al as a function of colors. Different symb ols indicate clu ster mem- 

bers selected using different membership criteria. The solid line indicates normal star colors from Kenvon & Hartmanr] dl995t> . 
while the dashed lines indicate the limits used to define the regions with IR excess stars {gray shade regions). Typical error bars 
in both directions are plotted in each panel. 



is found. Panel d shows the average r values of the three lines, 
where the dispersion is about 0.3. 

For the bulk of our stars, these values of r are consistent 
with zero or a small correction to the measured Li EW (consid- 
ering the dispersion as the typical error). In addition, the spread 
of the r values is comparable to the typical lithium EW error. 

We checked that r is not significantly correlated with either 
the IR excess indices or the FWZI of Ha. This suggests that our 
veiling determination is affected by large errors, and so would 
be for any correction to the Li EWs based on it. Such a cor- 
rection would have an error comparable to EW measurement 
errors themselves. We have applied this correction only for a 



few possible members with Li EW smaller than that in stars of 
similar temperature (see SectionHJl. 

4.3. Cluster membership criteria 

Our data allow us to assign the cluster membership to our tar- 
gets, based on the combination of RV measurements, X-ray de- 
tection and EW of the lithium line. In addition, the profile type 
of the Ha emission, when available, and optical-IR excesses, 
have been used as further indications of cluster membership. 

Of all 332 spectroscopically observed stars, we find 279 
single and 53 binary stars; of them 187 singles plus 52 binaries 
are X-ray detected, 190 single stars have RV consistent within 
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Fig. 13 Panels (a), (b) and (c) show the EW as a function of the temperature of the Vi, Nil and Fei lines used to estimate the 
veiling amount in our targets. The average values of the estimated veiling r as a function of temperature are shown in panel (d). 



3cr with the cluster mean RV, and 195 singles plus 53 binaries 
show a measurable lithium line. We consider the 158 X-ray de- 
tected single stars showing the lithium line, with a RV consis- 
tent with that of the cluster as certain cluster members. In addi- 
tion, we have 8 cases for which the Li EW is smaller than that 
in stars of similar temperature. We include them in the sample 
of certain cluster members; a discussion about these stars will 
be presented in Section[6] 

To the sample of cluster members we add 8 stars that are X- 
ray detected with a RV inconsistent with cluster membership, 
but Li EW larger than a temperature dependent threshold value 
(see Section [6] and Fig[9} and consistent with that of certain 
cluster members. We assume that these are binaries, members 
of the cluster, with the spectrum of only one of the two com- 
ponents being visible (SB1). We also consider as possible clus- 



ter members the 7 X-ray detected stars with RV inconsistent 
with cluster membership and Li EW smaller than the adopted 
threshold limit for the lithium EW; we labeled these stars with 
'M?'. Nevertheless, we cannot rule out that these are young 
field stars, unrelated to and slightly older than NGC 6530 stars. 

On the other hand, we have 10 additional single stars that 
are members based on the RV, with the lithium EW larger than 
the adopted threshold limit, but are not detected in X-ray. This 
is consistent with the expected incompleteness of our X-ray 
selected sample. All these stars have mass lower than 1.5 M 
(V > 16.3, Trff <4800K); we classify them as certain members 
('M') while we consider as possible cluster members ('M?') the 
3 stars that are X-ray undetected, with a RV consistent with the 
cluster members but with the lithium EW smaller than those of 



14 



L. Prisinzano et al.: VLT/Flames observations of the star forming region NGC 6530 



Table 5 Number of single X-ray detected and undetected stars 
(bottom Table) with different properties: RV members/non- 
members (Y/N) within 3cr, Li test certain members (Y), possi- 
ble members (Y?) and non-members (N) (see text). The num- 
ber of certain members are reported in boldface and those of 
possible members are in italic. 



Single X-ray detected 



\. Li 


Y 


Y? 


N 


Any 


Y 


158 


8 





166 


N 


8 


7 


6 


21 


Any 


166 


15 


6 


187 



Single X-ray undetected 



\. Li 
RV^^^ 


Y 


Y? 


N 


Any 


Y 


10 


3 


11 


24 


N 





1 


67 


68 


Any 


10 


4 


78 


92 



similar temperature stars; non members are indicated as 'NM'. 
Final membership results are reported in col. 20 of Table|3] 

The number of stars found with the adopted membership 
criteria are given in Table|5] We have a total sample of 184 
(158+8+8+10) single certain members (indicated in the table 
with boldface characters), where we also included the 8 candi- 
date depleted certain members (see Section[6j; by adding the 53 
binaries showing strong lithium 2 we have a total sample of 237 
certain members, plus 10 (7+3) possible single cluster mem- 
bers (indicated in the table with italic characters); thus, about 
74% of our initial sample are certain and possible cluster mem- 
bers. 

Our data also allow us to estimate the completeness of the 
X-ray data. By considering the sample of single stars, the X-ray 
completeness has been estimated as the ratio between the num- 
ber of X-ray detected stars, which are also members according 
to the lithium test and the total number of candidate members 
from Li, but regardless of X-rays. 

The fraction of X-ray detected stars in ACIS FOV and 
cluster locus, which were observed with GIRAFFE/FLAMES 
(setup HR15), is 0.61 (in the range V=14-18.2). The corre- 
sponding fraction for X-ray undetected stars (within the same 
magnitude and color range) is 0.25. Therefore, the inverse of 
these numbers (1.64 and 4.00, respectively) are the corrections 
for completeness, which we need to apply to our FLAMES 
samples to estimate the actual number of members (either X- 
ray detected or not) in the cluster locus down to V=18.2, which 
thus become (181*1.64)=296.66 (X-ray detected members) 
and ( 1 8 1 * 1 .64)+( 14*4.00)=352.66 (all members), respectively, 
if we include in these samples certain (M) and possible mem- 
bers (M?). Taking into account the contaminating members 
(6), the X-ray completeness is equal to 85% while, if we con- 
sider only certain members, the X-ray completeness is equal 
to 89%. The comparison of the two X-ray completeness val- 

2 all 53 binaries except one are X-ray detected 



ues allows us to estimate an error of the order of 4% in the X- 
ray completeness estimate. These values are consistent with the 
adopted X-ray incompletene ss correction applied to the Initial 
Mass Function in Prisinza no et al l il2005l) . 

5. Properties of IR excess stars: rotational 
velocities, accretion and binarity 

5.1. Introduction 

IR excesses are useful to discriminate PMS stars without disks 
(WTTS) from those with disks, which can also have accre- 
tion of material onto the star. On the other hand, T Tauri 
stars of similar spectral type show a trend with rotational ve- 
locity, in the sense that stars surrounded by accretion disks 
seem to be slow rotators. This has been explained in terms of 
disk locking jKoenigl[fl99TI) . but observational evidence for a 
correlation between rotat ion periods and accretion disk indi- 
cators has been debated JEdwards et all 19931: [Stassun et all 
ll999HHerbst et alll2000l) . In lRhode et all ll200ll) . the compar- 
ison of the vsin(i) distributions of CTTS and WTTS in the 
Orion Nebula Cluster, selected by means of / - K excesses 
and Ca II emission line strength, suggests, although the result 
is not statistically significant, that the CTTS sample contains a 
slightly larger fraction of slow rotators than the WTTS sample. 
On the contrary, using a different CTTS versus WTTS clas- 
sification based on the Ha line of the same sample of stars, 
Isicilia-Aguilar eta fl J2005I) found a significant difference be- 
tween the two distributions. More recently, using Spitzer data, 
a clear correlation betw een rotation a nd IR excesses has been 
found in Orion by Rebu ll et all (|2006). 

By using extinction-free indices Qabcd (with ABCD 
generic IR-optical magnitudes) as indicators of the presence 
of disks, we selected a sample of cluster members with IR ex- 
cesses, as described in Section lTTl On the other hand, we de- 
fined "No IR excess WTTS" as those members that do not show 
IR excesses in any of the considered indices Q. This allows 
us to select stars without strong disks and therefore, without 
strong signatures of accretion. In order to exclude high mass 
fast rotators, we considered only K and M type stars, with a 
temperature between 3700 K and 4500 K, estimated based on 
the V — I colors. 

5.2. Rotational velocities 

To quantify possible effects on the rotational velocities due to 
the circumstellar disks, we compared the distribution of the 
vsin(i) of stars with IR excesses with the distribution of the 
v sin(i) of the candidate WTTS stars, as shown in Fig.[H] K- 
S tests have been made to compare the v sin(i) distribution 
of the "No IR excess WTTS" and those of IR excesses stars. 
The resulting probabilities that the two distributions are dis- 
tinguishable (shown on the top of each panel) indicates that 
the distributions are statistically different at a significance level 
~ 98% for the Qvijk and Qbvjk excess stars and ~ 99% 
for the Qviu excess stars. This result allows us to support 
previous observational evidence that stars with disks (which 
can be CTTS stars) are slower rotators than "No IR excess 
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WTTS", as indicated in the disk locking pict ure and as al- 
ready found in the Orion Nebula Cluster jHerbst et all [2002; 
ISicilia-Aeuilar et al.Ll2005t lRebull et all bOOd) . 



Prob. 84.4% Prob. 97.8% 




50 100 ISO 200 00 100 150 200 

v sm(I) v s in(i) 

Fig. 14 Cumulative distributions of the v sin(i) values for the 
sample of the WTTS stars and for the samples of stars with 
IR excesses in the corresponding indices. The probability that 
the two distributions are significantly different is given in each 
panel. 




5 10 IS 20 25 5 10 15 20 25 

FWZKHn) [A] FWZl(H a ) [A] 



Fig. 15 Extinction-free indices of IR excesses vs. the Ha FWZI 
for the sample of stars for which the spectra in the Ha re- 
gion are available. Small dots and large bullets are WTTS and 
CTTS, respectively, defined based on the Ha line, X symbols are 
possible CTTS (see Section EOl and diamonds are stars with IR 
excesses in at least one of the four indices. The dashed lines in- 
dicate the limits used to define the regions with IR excess stars 
(gray shade regions). 

should be taken with caution. No further information can be 
deduced from Q indices for the possible CTTS. 



5.3. Ha 

For the subsample of members for which we also have Ha 
spectra, we plotted the Q indices versus the FWZI of Ha, as 
shown in Fig.^] where small dots and large bullets are WTTS 
and CTTS, respectively, defined based on the Ha line, while 
X symbols are possible CTTS (see Section l33t : diamonds are 
stars that show an IR excess at least in one of the four indices. 

29 out of the 3 1 stars classified as certain CTTS from Ha 
profile show IR excesses according to at least one index, and 
this is consistent with the presence of an accretion disk; on the 
other hand, we consider the WTTS stars with IR excesses as 
PMS stars with inert disks. The three stars classified as CTTS 
according to the Ha profile that do not show IR excesses are 
indicated with the 2MASS flag "c", which means that their 
photometry is biased by a nearby star that can contaminate 
the background estimation and therefore their IR magnitudes 



3 As already discussed, IR excesses are signatures of circumstellar 
material around the central star. However, not all stars with a disk 
show accretion phenomena and then our adopted indices cannot help 
us to distinguish CTTS with accretion disks from stars with a disk but 
without accretion. 



5.4. Binarity 

Fig-ED shows the optical-IR indices for binary stars. Dots 
are the SB2 stars, while triangles are the SB1 members (see 
Section l4~3l . Taking into account binaries with errors on each 
of the IR magnitudes smaller than 0.1 mag, we find that 10% 
(4/39) of them show Qjhhk excesses, 22% (9/41) show Qvijk 
and Qbvjk excesses, while only 6% (3/47) show Qviu excesses; 
all the 8 candidates SB1 are without IR excesses. The percent- 
age of SB2 stars with IR excesses in at least one of the four Q 
indices is 25% (10/40). 

On the other hand, among the single stars, 21% (29/137) 
show Qjhhk excesses, 51% (74/145) show Qvijk and Qbvjk ex- 
cesses, while 8% (13/160) show Qviu excesses. The percentage 
of single stars with IR excesses in at least one of the four Q in- 
dices is 57% (81/142). The result supports the hypothesis that 
compani ons do not contribute to the IR excess, as suggested 
byE ras-Velazauez & Ravi J2005I) . On the contrary, the signifi- 
cantly smaller fraction of binaries with disks (compared to sin- 
gle stars) indicates that binary stars in NGC 6530 have already 
undergone a significant disk evolution and inner disk clearing, 
probably due to gravitational effe cts of the binar y components 
on the circumstellar disks ( Armita ge et all 1 19991) . a theoretical 
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Fig. 16 Extinction-free indices of IR excesses for the binaries. SB2 stars are indicated with dots, while can didate SB1 s tars ar e 
indicated with triangles (compare with Fig. ll2t . The solid line indicates normal star colors from lKenvon & Hartmannl \ 19951) . 
The dashed lines indicate the limits used to define the regions with IR excess stars {gray shade regions). 



prediction for whi ch no clear observ ational evidence has been 
found until now dMonin et all [2006. and references therein). 

5.5. Summary 

We summarize in Table[6] the properties of member stars, ac- 
cording to the IR excess indices and the accretion indication 
from the Ho- line. In particular, in the upper Table, we give the 
number of X-ray detected and undetected members, the num- 
ber of members with RV at 3 cr consistent and inconsistent with 
cluster membership and the total number of members, for the 
stars with IR excesses in the specified index. The numbers cor- 
responding to the "Q xxxx " line indicate the total members hav- 
ing IR excess in at least one of the four indices, while the num- 



bers corresponding to the line "All" indicate the total numbers 
of members, regardless of their photometry. 

Analogously, in the lower Table, we give the same numbers 
with respect to the sample of members with Ha observations. 

The percentage of members recovered by X-ray detection 
is 95% (174/184), while the members showing IR excesses are 
44% (81/184). Among X-ray missed members, 7 out of 10 are 
recovered using only optical-IR photometry; this means that 
the percentage of members not recovered by X-ray detections 
and/or IR excess diagnostics is equal to about 2% (3/184). 

If we consider the sample of cluster members with Ha ob- 
servations, we find that approximately 40% of the stars show 
accreting disks (profile types PCI, PCIII, YY Ori in Table Q). 
This is a percentage similar to that found in the Orion Nebula 
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Table 6 Top: Number of X-ray detected and undetected mem- 
bers, with RV at 3 <x consistent and unconsistent with cluster 
membership and total number of members, for the stars with IR 
excesses in the specified index. The numbers corresponding to 
the "Q AAAA " line, indicate the total members having IR excess in 
at least one of the four indices, while, the numbers correspond- 
ing to the line "All" indicate the total of members, regardless 
their photometry. Bottom: same numbers with respect to the 
sample of members with H„ observations. 



Total single member sample 





X-ray 


No X-ray 


RV 


NoRV 


Tot. 


QjHHK 


25 


4 


29 





29 


QviJK 


67 


7 


74 





74 


Qvm 


11 


2 


13 





13 


Qbvjk 


68 


6 


74 





74 


Qxxxx 


74 


7 


81 





81 


All 


174 


10 


176 


8 


184 


Single member sample with H„ obs. 




X-ray 


No X-ray 


RV 


NoRV 


Tot. 


CTTS 


29 


2 


31 





31 


All 


67 


4 


69 


2 


71 



Cluster bv lSicilia-Aguilar et alJ (Eo05) with similar indicators. 



6. Lithium abundances and depleted members 

As mentioned in Section l4~3l we consider as cluster members 
all objects with lithium EW larger than the adopted threshold 
(see dashed line in Fig.|9jl, having a RV consistent with the 
cluster membership and/or an X-ray counterpart. The adopted 
membership criteria always require the presence of the lithium 
line, because this condition has to be always true for young 
stars, such as those of NGC 6530. In the sample of certain clus- 
ter members we also included X-ray detected stars, with RV 
consistent with the cluster membership and EW of the lithium 
measurable but smaller than the EW threshold. 

While the lithium EWs strongly depend on temperature, the 
lithium abundance is expected to be constant for very young 
stars such as our targets, which should not have yet reached 
interior temperatures for lithium burning. 

However, accurate lithium abundances cannot be derived 
from our data due to the large uncertainty affecting both our 
temperature estimates and lithium EW measurements. The 
large spread in temperature is due to the large errors mainly 
introduced by unknown individual reddening values, while the 
spread in the EWs is mainly due to the continuum normaliza- 
tion and to the veiling, which has been estimated to be of the 
same order as the random errors affecting the EW measures. 

In addition, the growth curves given in the literature to de- 
rive lithium abundances from E Ws are avai l able o nly for stars 
with 4000 <T< 6500 K fcoderblom et all Il993h . while our 
sample contains stars with temperatures down to T~3500K, 
where the physics of stellar atmospheres is not well known due 



to the presence of molecular lines such as TiO and H2O. A 
further uncertainty is introduced by non-lo cal thermodynamic 
equilibrium (NLTE) JCarlsson et all 1 1994 . 

Using our data, the Soderblom et al. I ill 9931) growth curves 
and the NLTE lithium correction of ICarlsson et al.l (11994ft . we 
find that, for the sample of stars classified as members and 
T > 4000 K, the median value of the lithium abundance is 
A(Li)=3.0 dex 4 , with a standard deviation of 0.6 dex. While the 
median value is consistent with the expected cosmic lithium 
abundance, the spread is quite large and is mainly due to the 
few objects that have a lithium EW smaller than those of simi- 
lar temperature stars. Based on the EW values of the bulk of 
our data, and taking the EW of the Pleiades as a reference 
for older stars, we consider the stars with EWs larger than 
the dashed line in Fig. [9] and temperature smaller than about 
5000 K as lithium-undepleted cluster members. We also con- 
sider as lithium-undepleted members the few objects with tem- 
perature higher than 5000 K. 

Particular attention has been devoted to the X-ray detected 
objects, with a RV consistent with cluster membership and EW 
smaller than the limit indicated by the dashed line, as shown in 
Fig.|9l Tabled lists the values obtained in this paper for these 
stars; the estimated veiling corrections (r) are also reported for 
the stars in which the three lines adopted to estimate the veiling 
are visible. Using these values to correct our measured EW of 
the lithium line, we find that the corrected EWs are consistent 
with the values of the Li EW of similar temperature members 
for the stars with spectrum numbers Sp 1128, 3042 and 3099, 
while, for the star Sp 2079, it remains smaller than the adopted 
threshold. Note that the r values are among the highest values 
in Fig. The age estimated for the star Sp 2079 (15Myr) is 
consistent with a significant lithium depletion, as predicted by 
theoretical models (see Fig. ll7l >. However, this value can be af- 
fected by a systematic error due to the veiling, which can mod- 
ify the photospheric colors of the stars. In addition, its temper- 
ature, among the highest in our sample, is affected by a large 
error due to the unknown reddening. We do not exclude the 
hypothesis that this object is a young main sequence field star. 

For the remaining 4 stars, only the lithium line is clearly 
visible and this suggest that they are highly veiled objects. This 
hypothesis can be confirmed for the stars 2106 and 2110, for 
which evidence of accretion has been found from the Ha line. 
The presence of circumstellar disks, deduced by IR excesses, 
can be considered as a further indication of possible accretion 
phenomena which can cause significant veiling. 

The position of all our sample stars in the observed V vs. 
V - I color-magnitude diagram is shown in Fig.[l7] where 
the de reddened theoretical isochrones and tracks of lSiess et alJ 
(2000) are drawn. The position of the 8 stars with low Li EWs 
and the regions with expected lithium depletion are also indi- 
cated; portions of the tracks with undepleted Li are drawn as 
dotted lines, those with predicted Li depletion up to a factor of 
10 are drawn as dotted-dashed lines and those with predicted 
Li depletion more than a factor of 10 are drawn as solid lines. 

Temperature estimates from spectral types, allowing a con- 
sistent comparison with theoretical models, will be helpful to 

4 A(Li)=log(n L , /«„)+! 2 lAnders & GrevesseU 19891) 
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establish whether these stars are actually lithium depleted or if 
the EW of these stars are smaller than those of similar temper- 
ature stars because of only veiling effects, or whether they are 
young field stars. 




1.0 1.5 2.0 2.5 3.0 3.5 

V-I 



Fig. 17 V vs. V-I color-magnitude diagram of candidate mem- 
bers of NGC 65 30. Theoretical isochrones and tracks of 
ISiess et alJ §000) are drawn and different levels of predicted 
Li depletion are also indicated; portions of the tracks with un- 
depleted Li are drawn as dotted lines, those with predicted Li 
depletion up to a factor of 10 are drawn as dotted-dashed lines 
and those with predicted Li depletion more than a factor of 10 
are drawn as solid lines. The 8 stars with evidence of Li deple- 
tion and/or veiling effects are also indicated. 



7. Summary and conclusions 

We used spectra of intermediate resolution (R~19 300), taken 
with the GIRAFFE-FLAMES multi-object spectrograph of the 
ESO-VLT Kueyen Telescope, to analyze a sample of 332 can- 
didate members of the young cluster NGC 6530. The targets 
were selected in the region of PMS stars of the V vs. V-I 
color-magnitude diagra m, traced by a large number of X-ray 
detected stars (Prisinzan o et allEoO SL 72% of the selected tar- 
gets are X-ray active stars. 

The spectra of these stars allowed us to measure both ra- 
dial and rotational velocities using the technique of Cross- 
Correlation. This technique enabled us also to identify 53 ob- 
jects as binaries. In addition we measured the EW of the lithium 
line, used as a membership criterion. Ha spectra, available for 
a subsample of 1 15 stars, allowed us also to select WTTS and 
CTTS. 

The measures obtained from our spectra were used to se- 
lect a sample of certain and possible members, based on the 



combination of the RV measures, the Li EWs and the X-ray 
detection. We find 176 certain single members that, added to 
the 53 binaries, gives us a total sample of 237 certain members 
plus 10 possible members. 

We used extinction-free indices to select stars with IR ex- 
cesses, indicators of the presence of circumstellar disks. A 
statistically significant difference was found between the ro- 
tational velocity distributions of K and M type stars with IR 
excesses with respect to that of the same spectral type stars 
but without IR excesses. In addition we find that all stars, ex- 
cept three, classified as CTTS based on the Ha line have IR 
excesses and this confirms that most of the IR excesses stars 
have accretion phenomena. These results are consistent with 
the disk locking picture suggesting a rotational velocity reduc- 
tion in st ars with disks, as already found in the Orion N ebula 
Cluster JSicilia-Aguilar et all 120051: iRebull et all 120061 and 
reference therein). 

We also find that the fraction of binaries with disks is sig- 
nificantly smaller than that of single IR excess stars and this 
suggest that gravitational effects of the binary components in 
NGC 6530 caused a significant disk evolution and disk clear- 
ing. 

A fraction equal to about 44% (81/184) of the stars iden- 
tified as cluster members show IR excesses and, among the 
X-ray undetected members, about 70% (7/10) have IR excess, 
recovered using only optical-IR photometry. This confirms the 
importance of using both X-ray data and optical-IR data as di- 
agnostics of membership. 

The EWs of the lithium line of the stars classified as clus- 
ter members are significantly larger than those of older stars 
(Pleiades like) of similar temperature and this is agreement 
with the young ages of NGC 6530. We computed the lithium 
abundances of cluster members; while the median value of 
these values, equal to A(Li)=3.0 dex, is in agreement with the 
predicted cosmic abundance, the rms, equal to 0.6 dex, is quite 
large; this is due to the uncertainty in the computed abundance 
values due to the uncertainties affecting the temperature esti- 
mates, the lithium EW measures (including the veiling correc- 
tion), the adopted growth curves and the non-local thermody- 
namic equilibrium. 

We find 8 possible cluster members that are X-ray detected 
stars and with RV consistent with cluster membership, with Li 
EW smaller than the adopted threshold. For 7 of them, we find 
evidence of strong veiling, which can explain their small Li 
EW values, while for one we find that the veiling corrected EW 
is however smaller than those of stars of similar temperature. 
The age estimated for this star is consistent with a significant 
lithium depletion, but we do not exclude that possible system- 
atic errors due to the veiling can affect the determination of 
stellar parameters (including age). 
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Table 7 Properties of stars with evidence of Li depletion and/or veiling effects. 
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Age 


M 


vsin(i) 


G~ vsin(i) 


T e ff 


Li i 






X-ray 


RV 


Q 


Q 


Q 


Q 


H a 


H a 


Type 




(Myr) 


(M a ) 


[km/s] 


[km/s] 


K 


EW[mA] 




FWZI[A] 






JHHK 


VIJK 


VIIJ 


BV.IK 


CTTS 


PT 




1128 


1.06 


0.61 


21.9 


2.0 


3902 


420. 


0.3 


_ 


Y 


Y 


Y 


Y 


Y 


Y 


_ 


_ 


M 


2079 


15.03 


1.11 


10.8 


0.5 


4648 


176. 


0.5 


2.2 


Y 


Y 






N 




WTTS 


AE 


M 


2106 


0.44 


0.81 


30.2 


4.6 


4108 


268. 


1.9 


15.6 


Y 


Y 


Y 


Y 


N 


Y 


CTTS 


PCI 


M 


2110 


1.35 


1.06 


59.6 


7.7 


4319 


307. 


NaN 


13.2 


Y 


Y 


Y 


Y 


Y 


Y 


CTTS 


PCIII 


M 


3013 


5.42 


0.94 


< 10.0 




4190 


220. 


NaN 




Y 


Y 


N 


Y 


Y 


Y 






M 


3042 


1.12 


0.79 


< 10.0 




4063 


366. 


0.3 




Y 


Y 


Y 


Y 


Y 


Y 






M 


3099 


9.48 


1.00 


35.2 


5.1 


4272 


345. 


0.7 




Y 


Y 






N 








M 


3109 


5.47 


0.85 


< 10.0 




4104 


322. 


1.9 




Y 


Y 


N 


Y 


N 


Y 






M 



Typical error on the veiling correction is 0.3 
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